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An  attractive  approach  to  short-range  forecasting  is  to  determine  cloud 
motion  from  a  sequence  of  satellite  images  and  extrapolate  the  patterns  and 
associated  weather  into  the  future.  Objective  motion  vector  techniques  are 
available  and  the  forecast  procedure  can  be  accomplished  by  computer.  This 
approach  is  being  evaluated  at  AFGL  and  this  report  presents  results  of 
testing  motion  vector  techniques.  Tracking  and  covariance  techniques  were 
,  compared  along  with  winds  aloft  and  persistence  (no  motion)  as  controls.  A 
f  covariance  technique  had  top  score,  but  only  slightly  better  than  persistence 
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LOCAL  FORECASTING  THROUGH  EXTRAPOLATION  OF 
COES  1MACERY  PATTERNS 


II.  Stuart  Muench 

Air  Force  Geophysics  Laboratory 
Bedford,  Massachusetts 


1 .  INTRODUCTION 

The  development  of  the  geo¬ 
synchronous  weather  satellite  (GOES)  has  pre¬ 
sented  forecasters  with  the  opportunity  to 
significantly  advance  the  art  of  short-range 
forecasting.  Where  the  forecaster  formerly 
had  hourly  surface  observations  from  stations 
100  km  apart,  he  now  can  have  digital  satellite 
observations  every  half-hour  at  locations  1-km 
apart  for  video  and  8-km  apart  fob-  IR.  In 
fact,  the  previous  famine  for  mesoscale  data 
lias  now  become  a  feast.  A  single  satellite 
transmits  more  numbers  in  a  few  seconds  than  do 
all  the  weather  observers  in  the  world  in  24 
hours.  The  problem  for  the  forecaster  Is  now 
how  to  digest  all  this  information  and  still 
meet  the  deadlines  of  short-range  forecasts 

An  attractive  way  to  assist  the 
forecaster  is  to  develop  an  automated  (or  semi- 
automated)  forecast  procedure  which  produces 
guidance  forecasts.  An  elementary  procedure, 
based  on  extrapolation,  is  shown  below  in 
Figure  1. 
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Fig.  It  Illustration  ol  short-range  fore¬ 
casting,  using  GOES  Imagery  data. 


At  the  top  of  Figure  1  are  two 
overlaid  satellite  Images,  the  earlier  image 
shown  by  dotted  lines.  One  can  use  objective 
or  subjective  techniques  to  obtain  a  motion 
vector,  as  indicated  by  the  short  arrow 
pointing  towards  station  XXX.  To  forecast 
for  XXX,  one  reverses  the  vector  and  looks 
"upstream"  Increments  of  distance  in  pro¬ 
portion  to  the  increments  of  time  that  are  of 
interest.  At  each  Increment  one  extracts  the 
satellite  observed  reflectivity  from  the  most 
recent  image,  and  the  result  is  the  forecast 
of  reflectivity  in  Figure  IB.  Finally  one 
uses  statistically  derived  algorithms  to  con¬ 
vert  reflectivity  to  cloud  cover  and  pro¬ 
bability  of  precipitation  as  shown  in  Figure  1C. 

This  approach  has  been  used 
with  digital  radar  data,  both  by  the  National 
Weather  Service  and  McClll  University,  and 
should  be  applicable  to  GOES  Imagery  data. 

A  more  sophisticated  version  would  Include 
Infrared  as  well  as  video  reflectivity,  and 
more  forecast  parameters.  However,  we  should 
first  test  the  concept  with  a  simple  model  to 
determine  if  changes  In  reflectivity  (hence 
weather)  can  be  adequately  forecast  and 
determine  where  sophistication  is  needed  to 
improve  performance. 

2.  FIRST  TEST  OF  MOTION  VECTOR 

TECHNIQUES 

Several  objective  techniques 
for  determining  awtlon  vectors  from  successive 
satellite  Images  have  been  reported  In  the 
literature.  In  general,  these  were  developed 
to  be  operated  Interactively  with  a  computer 
to  measure  winds  aloft  in  regions  of  sparce 
data.  However,  the  cloud  tracking  procedures 
of  SRI  (Endlich  et  al-1971,  Wolf  at  al-1977) 
and  the  cross-covariance  technique  (Leese  and 
Novak-1971)  appear  capable  of  operating  in  an 
automatic  mode.  A  copy  of  the  SRI-1977  program 
was  obtained  from  the  Naval  Environmental  Pre¬ 
diction  Research  Focllity  and  computer  programs 
were  constructed  based  on  journal  articles  for 
the  SRI-1971  and  the  cross-covariance  technique 
(using  fast-Fourter-Transforms).  While  these 
programs  were  being  prepared,  a  binary  co¬ 
variance  technique  was  developed,  which  reduces 
the  GOES  (mage  to  binary  values  (one  or  aero) 
and  computes  covariance  using  logical  arith¬ 
metic,  bO  grldpolnts  at  a  time. 


The  following  are  characterletlce 
of  Che  four  motion  vector  techniques  (more  de- 
talla  can  be  found  In  Muench  and  Hawkins,  1979)i 

1)  SRI-1971 

30x30  arrays,  eelect  100  brlghteet  point*. 
Identify  and  locate  cloud  centere. 

Hatch  pa  Ire,  compute  mean  displacement. 

2)  SRI-1977 

70*70  array*,  eelect  brlghtnee*  6  percent, 
Locate,  Identify  and  characterlae  cel  la 
Hatch  pair*,  compute  mean  dleplacement. 

3)  I Ft  Croae-covar lance 

32x32  arrays,  convert  to  frequency  domain. 
Compute  covariance  matrix  of  trial 
dleplacement*.  Find  dlaplacement  with 
maximum  covariance. 

A)  Binary  Croae-covarlanc* 

60x60  Initial  and  80x80  final  arrayt. 
Convert  to  binary.  Sum  matching  l'e 
and  O'*  for  trial  dleplacement a. 

Find  dleplacement  with  maximum  match**. 

The  flrat  text  uaed  a  240x240  array 
of  2-km  COES  video  data  aa  In  initial  array. 

For  the  final  array,  the  Initial  Imagery  wax 
uaed  but  artificially  dlaplaced  to  almulate  a 
motion  of  13  mp*  over  $  hour.  Each  motion 
vector  program  extracted  data  value*  ao  aa  to 
cover  the  array*  (e.g.  every  7th  point,  every 
3rd  point).  Six  Image*  were  choxen,  repre¬ 
senting  a  variety  of  synoptic  conditions, and 
for  each  image,  3  different  displacements  were 
made.  In  all,  18  motion  vector*  were  computed 
for  each  technique.  Some  preliminary  testa 
Indicated  that  the  SRI  programs  were  having 
difficulty  with  fine  structures,  and  so  the 
240x240  arrays  were  given  an  8xl6-polnt  smooth¬ 
ing.  The  results  of  the  test  are  shown  below 
In  Table  1. 

Table  1.  Test  results  using  artificial 
displacements.  Vector  error 
Is  rms  resultant  error  divided 
by  artificial  displacement. 


rms  col. 

error 

rms  row 

error 

Vector 

error 

SRI-1971 

♦8.3 

♦3.9 

♦59% 

SRi-1977 

♦6.0 

♦  3.7 

♦49% 

fFt  Cross 
covariance 

♦4.9 

♦  3.8 

♦42% 

Binary  Cross- 
Covariance 

♦0.3 

♦  1,3 

♦  10% 

In  this  first  test,  the  binary 
cross-covariance  technique  was  clearly  super¬ 
ior  to  the  other  three.  An  inspection  of 
individual  cates  showed  that  all  techniques 
worked  well,  vector  errors  I  St  or  less,  when 
the  brightest  portion  of  the  clouds  was  in 
the  center  of  the  array.  For  a  variety  of 
reasons,  the  first  three  techniques  produced 
unreliable  results  when  the  brightest 
portion  was  on  a  boundary. 


3.  SECOND  TEST  OF  MOTION  VECTOR 

TECHNIQUES 

The  first  test  was  well  con¬ 
trolled  in  that  the  true  motions  were  pre¬ 
determined.  And,  the  test  was  revealing 
since  the  boundary  problems  were  uncovered. 
However,  the  test  was  not  completely  realistic 
aa  the  consecutive  Images  did  not  Include  the 
effects  of  cloud  development  and  decay,  as  well 
as  motion.  Thus,  a  second  test  was  run  to  check 
the  performance  of  the  techniques  using  series 
of  consecutive  half-hourly  video  Images.  ■  The 
object  of  the  test  was  to  have  each  of  the  mot¬ 
ion  vector  techniques  determine  motion*  and 
make  forecasts  of  reflectivity  through  simple 
extrapolation  of  pattern  motion.  Verifications 
of  forecasts  would  Indicate  which  techniques 
produce  the  most  reliable  result*  and  com¬ 
parisons  with  a  control  such  as  persistence 
would  indicate  to  what  extent  simple  extra¬ 
polation  it  a  useful  short-range  forecast  pro¬ 
cedure  . 

Twelve  cates  were  selected, 
representing  widely  varying  synoptic  conditions, 
including  cold  fronts,  warm  fronts,  advancing 
and  rttreatlng  cyclonic  storms.  Each  cate 
consisted  of  tlx  consecutive  satellite  video 
images  at  normally  half-hourly  Intervals  (In 
four  cate*  there  was  a  one-hour  gap  between  two 
of  the  Images).  The  region  for  which  2-km 
video  dtta  are  archived  Is  shown  below  In 
Figure  2. 


Figure  2:  Satellite  data  archive  area,  location* 
of  motion  vector  points,  and  grid 
of  points  used  for  forecast  and 
verification. 

Four  techniques  were  added  to 
those  used  in  tha  first  test.  A  variation  was 
made  on  the  fFt  cross-covariance  to  have  a 
larger  48x48  initial  array  In  hopes  of  att- 
tvlatlng  the  boundary  problem,  through  risking 
aliasing  problems.  For  controls,  two  radio¬ 
sonde  wlndt-aloft  were  included,  the  1200  I'T 
700mb  level  wind  and  4  the  SOOmb  level  wind. 

These  winds  tr*  sometimes  used  to  "steer"  radar 
pattern*  and  developing  storms.  The  last  con¬ 
trol  added  was  "persistence"  or  no  local  change 
with  time,  a  control  that  It  often  difficult  to 
beat  in  short  rang*  forecasts. 


2 


..  The  five  notion  vector  programs 
using  satellite  data  computed  motion  vectors 

at  each  of  the  four  locations  shown  In  Figure  2, 
for  each  consecutive  pair  of  images  in  each  of 
the  12  cases*  A  25-pt  array*  with  6A-km 
spacing,  was  set  up  over  each  location  to 
serve  as  the  forecast  and  verification  points. 


Table  2.  Percent  Correct  and  Skill 

Scores  of  Forecasts  from  Motion 
Vector  Techniques.  Twelve 
cases,  all  time  intervals 
and  all  thresholds  combined. 


In  order  to  assure  correct  posi¬ 
tioning  of  the  satellite  data,  a  two  stage 
"navigation"  procedure  was  used.  The  basic 
satellite  attitude  and  orbital  parameters 
are  refined  dally,  using  NESS  bulletins  and 
minimising  differences  between  predicted  and 
observed  locations  of  prominent  geographic 
features  (usually  in  Mexico,  Venesuela  or 
Bolivia).  This  navigation  procedure  is  carried 
out  on  the  AFGL  MclDAS  system(computer  inter¬ 
active  graphics),  and  overall  gives  positioning 
accuracy  of  about  ^20km  in  the  NE  United 
States.  A  "fine  tuning”  of  each  archived  image 
is  performed  using  identifiable  landmarks,  with¬ 
in  the  region  shown  In  Figure  2,  which  reduces 
the  positioning  error  to  about  +5  km. 

Varying  sun  angle  causes  pro¬ 
blems  when  one  uses  satellite  video  data 
quantitatively,  and  measures  were  taken  to 
correct  for  sun  angle.  Within  a  case,  each 
image  was  normalised  to  noon  (EST)  of  Chat  a 
day,  through  use  of  a  factor  determined  by  com¬ 
paring  the  histograms  of  reflectivity.  Four 
reflectivity  thresholds  were  selected  to 
verify  forecasts,  0.30,  0.55,  0.75,  and  0.90, 
normalised  to  overhead  sun.  For  each  case, 
the  thresholds  were  adjusted  to  the  noontime 
sun,  using  an  anisotropic  scattering  routine 
described  in  Muench  and  Keegan  (1979). 

The  forecasts  and  verification 
were  made  by  reversing  the  specified  motion 
vector  and  computing  the  "upstream"  locations 
at  intervals  equivalent  to  half-hour  linear 
motion,  out  to  24  hours.  The  forecast  reflect¬ 
ivity  consisted  of  an  average  of  the  four 
nearest  2 -km  resolution  satellite  values 
(normalised  to  local  noon).  With  respect 
to  each  of  the  thresholds,  forecast  and 
verification  reflectivity  were  either  below 
or  else  equal  to  or  above  the  threshold.  Scores 
were  kept  of  the  numbers  of  times  both  fore¬ 
cast  and  verification  were  In  the  tame  category 
and  numbers  of  times  they  were  not.  The  scores 
were  summed  over  all  cases  and  final  results 
presented  in  terms  of  percent  correct.  In 
addition,  skill  score  relattve  to  persistence, 
SSp,  was  computed  byi 

SS  «(F-P)/(1.-P) 

P 

wehre  F  Is  percent  correct  for  the  motion  vector 
technique  and  P  is  the  percent  correct  for  per¬ 
sistence.  Results  summed  over  all  forecast  time 
Intervals  and  all  thresholds  are  shown  below 
In  Table  2  (More  completed  results  can  be 
found  in  Muench-1979). 


Technique  Percent  Correct  Skill  Scores 

Satellite-based  motion  vectors 


SR1-197KAFCL) 

0.880 

-0.16 

SRI-1977 

0.892 

-0.05 

fFt(l)  Covariance 

0.897 

0.00 

f Ft (2)  Covariance 

0.900 

0.03 

Binary  Covariance 

0.903 

0.06 

Controls 

700 -mb 

0.900 

0.03 

4  500  mb  wind 

0.898 

0.01 

Persistence 

0.897 

0.000 

Amont  the  satellite-based  tech¬ 
niques,  the  relative  standings  are  exactly  as 
they  were  in  the  first  test,  with  the  best  score 
for  the  binary  covariance,  followed  by  the  fFt- 
covarlance  techniques,  then  the  SRI-1977  and 
finally  the  SRI-1971.  In  Figures  3  and  A 
(following  page)  the  relative  positions  are 
essentially  the  same  for  all  time  intervals  and 
all  thresholds,  so  even  though  the  differences 
are  small,  the  results  appear  stable.  In 
the  lower  part  of  Table  2,  one  notes  that  the 
two  upper-wind  controls  scored  nearly  as  well  as 
the  covariance  techniques,  raising  the  question 
whether  the  effort  of  computing  covariance  motion 
vectors  can  be  justified  in  regions  where  radio¬ 
sonde  data  are  plentiful. 

The  most  Important  point  In 
Table  2  and  Figures  3  and  A  is  that  none  of  the 
techniques  Is  much  better  than  simple  persis¬ 
tence.  This  becomes  readily  apparent  by 
examining  the  skill  scores  in  Table  2,  which 
Indicate  some  are  actually  negative.  While 
the  percent  correct  for  persistence  is  quite 
high,  leaving  little  room  for  improvement, 
never-the-less,  the  best  of  the  motion  vector 
techniques  correctly  forecasted  less  than  one 
of  every  sixteen  threshold  crossings  that 
occurred. 
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The  firtt  thought  upon  sating 
tho  lull  tklll  ocoroo  was  that  thoro  warn  pro¬ 
gramming  or  navigation  errors •  However,  a 
thorough  chock  ahowad  that  thla  was  not  tha 

caaa. 

Next,  thara  was  eoncarn  over  tha 
uaa  of  4x4km  (2x2  milt)  avarag*  reflectivity 
in  tha  foracatta  and  verification.  Moat  cloud 
calla  at  email  at  4km  could  not  ba  expected  to 
have  a  ltfe-apan  more  than  hour,  and  thoae 
that  did  would  require  very  praclte  determination 
of  motion  to  ba  correctly  forecaet(  particularly 
when  apeeda  were  greater  than  lOnpt.  So,  tha 
foracaet-verlflcatlon  program  waa  re-run,  only 
a  23x23km  (13x13  mile)  emoothing  waa  applied 
to  the  reflectivity  data.  However,  the  recultlng 
acoree  were  again  nearly  the  aamc  aa  peralatence. 
Figure  3  ahowa  the  variation  of  percent  correct 
with  time  for  both  peralatence  and  the  binary 
croas-covarlance  technique,  with  the  lower 
curves  for  4x4km  run  and  the  upper  curvea 
for  the  23x23km  run.  The  principal  effect  of 
the  smoothing  was  to  Increase  the  percent  correct 
for  persistence  at  all  time  Intervals,  with  a 
nearly  identical  increase  for  the  covariance 
technique.  Apparently,  the  larger  disturbances 
remaining  after  the  23km  smoothing  were  just 
as  difficult  to  forecast  as  the  smaller  ones 
filtered  out. 


'  Figure  J.  Percent  correct  tor  forecasts 

by  time  intervals  (4  thresholds 
combined). 
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Figure  4.  Percent  correct  for  forecaats, 
by  threshold  <3  time  intervals 
combined). 
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Figure  3.  Percent  correct  versus  time,  for 
binary  cross-covariance  and  per¬ 
sistence,  for  4x4km  average  data 
and  23x23km  average  data. 
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DIAGNOSIS  OF  PROBLEMS 
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An  obvious  conclusion  from  the 
second  test  is  that  the  simple  extrapolation 
forecast  scheme  lllustrared  in  Figure  1  will 
not  produce  forecasts  substantially  better 
than  persistence.  While  we  might  consider 
trying  to  Improve  the  motion  vector  techniques, 
at  this  point  It  would  be  prudent  to  perform 
a  more  complete  diagnosis  to  see  If  perhaps  the 
scheme  is  too  simple  and  that  more  sophistica¬ 
tion  is  required. 


An  examination  was  made  of  a 
case  on  IS  Nov,  1977,  A  weak  cold  front  had 
passed  through  the  northeast  the  previous  day, 
leaving  an  eastwest  quasl-stat ionary  front 
across  the  Virginias.  While  there  was  still 
cold  air  advection  in  the  lower  levels  over 
New  England,  a  weak  upper  level  trough  was 
moving  rapidly  eastward  from  the  Mississippi 
Valley,  along  with  a  slowly  developing  wave 
on  the  front. 


Normalised  reflectivity  values 
were  extracted  for  31x41  arrays  of  10km 
average  data  over  each  of  the  stations  shown 
In  Figure  2.  Similar  arrays  were  also  obtained 
lor  lnirared  "temperature"  and  reflectivity 
standard  deviation.  Surface  observations 
and  radiosonde  ascents  were  analysed  to  deter¬ 
mine  the  three-dimensional  cloud  structure  and 
its  variations  with  time. 

An  example  of  an  analysis  of 
the  reflectivity  field  around  station  ALB  is 
shown  below  in  Figure  6.  In  the  lower  right, 
the  skies  were  essentially  clear,  while  to  the 
northwest  clouds  prevailed.  Over  the  hills  of 
Vermont  and  New  Hampshire  were  heavy  strato- 
cumul 1  that  were  stationary,  such  as  the 
stlpled  areas  marked  "A”.  However,  a  clear 
area  in  northeast  New  York  state,  marked  ”B" 
was  spreading  steadily  south  and  east  during 
the  day  and  eventually  the  clouds  over  the 
mountains  disappeared.  Obviously,  no  single 
motion  vector  could  correctly  forecast  the 
features  A  and  B. 


Figure  hi  Analysis  of  reflectivity  in 

vicinity  of  station  ALB,  at  1700  I'T, 


To  the  west  of  station  ALB, 
another  complication  was  found.  The  heavy 
stipled  areas  narked  "C”  were  the  most  Intense 
portion  of  a  stationary  pattern  that  extended 
southeast  from  Lake  Ontario.  These  were 
basically  low  clouds  with  considerable  vertical 
development,  including  rain  showers.  Above, 
there  was  a  layer  of  cirrus  that  had  spread 
eastward  to  near  ALB,  and  lnbedded  in  the 
cirrus  was  a  dense  area  marker  "D"  that 
was  moving  east  at  30  to  33  taps.  Again, 
no  single  motion  vector  could  forecast  these 
Independent  motions. 

Figure  7,  below,  shows  another 
analysis  of  the  reflectivity,  near  station 
WAL  in  the  Del-Mar-Va  peninsula.  In  this 
region  there  was  an  extensive  area  of  middle 
and  high  clouds  that  was  spreading  eastward, 
and  rain  was  associated  with  che  stipled  area 
clouds.  A  few  holes  in  the  clouds  were  present 
in  the  upper  right  (northeast)  portion.  The 
reflectivity  pattern  was  quite  complex  in 
form,  and  as  might  be  expected,  the  forms 
changes  from  image  to  image.  Part  of  the 
change  was  due  to  a  systematic  weakening  of 
the  most  Intense  areas  as  they  moved  east- 
northeastward.  Other  changes  have  not  yet 
been  diagnosed.  Such  changes  would  lower 
the  scores  of  extrapolation  based  forecast 
techniques.  Also,  within  this  area,  the 
patterns  in  the  upper  left  were  moving 
noticeably  faster  than  those  in  the  lower 
right,  and  a  single  motion  vector  was  not 
adequate. 


Figure  7i  Analysis  of  reflectivity  in  vicinity 
of  station  WAL,  at  1700  UT,  15  Nov 
1978. 

At  more  cases  are  analysed, 
undoubtedly  other  complicating  factors  will 
be  found.  However,  at  this  point  it  it  clear 
that  a  successful  forecast  technique  will  have 
to  Include  a  means  to  cope  with  patterns  with 
different  motion,  close  to  each  other  or  even 
superimposed.  A  partial  solution  is  to  per¬ 
form  a  vertical  separation  of  the  cloud  layers, 
using  all  available  data,  and  forecast  each 
layer  separately. 


W«  thou Id  keep  In  alnd  th*t 
th«  satellite  imagery  fro*  clouds  i«  airily 
a  manifestation  of  vortical  air  notion  in  a 
Mist  environment.  Tha  cloud  partlclaa  grow 
and  fall  aufflciantly  faat  that  images  an 
hour  or  more  apart  ara  not  viawtng  tha  seme 
partlclaa.  In  aaaanca.  whan  wa  foracaat 
tha  futura  cloud  pattarna,  wa  ara  foracaating 
tha  varttcal  notions .  So,  ona  would  do  wall 
to  kaap  tha  dynaaiica  and  phyaica  in  nind 
whan  davaloplng  satellite-based  foracaat 
techniques. 

3.  SUMMARY  AMD  CONCLUSIONS 

To  utlliaa  gaoaynchronoua 
aatalllto  information  in  ahort-ranga  fora¬ 
caating,  a  almplc  axtrapolation  procadura 
waa  propoaad.  Succaaaiva  imagery  data  would 
ba  uaad  to  conputa  a  motion  vactor.  Next, 
the  motion  vactor  would  ba  used  to  pradlct  tha 
time  hiatory  of  raflactivlty  at  a  point. 
Finally,  algorithms  would  translate  the 
reflectivity  forecast  to  common  weather 
parameters. 

Test  were  run  to  find  the 
most  suitable  means  to  compute  motion  vectors. 
In  the  first  test,  using  specified  motions, 

Che  techniques  described  in  the  literature 
worked  well  only  when  the  brightest  portion 
of  the  pattern  was  in  the  center  of  the 
working  array.  A  binary-cross-covariance 
technique  worked  well  for  all  conditions.  In 
a  second  test,  making  reflectivity  forecasts, 
again  the  binary  cross-covariance  scored  beat, 
but  the  differences  were  small.  The  scores 
of  the  upper-level  winds  were  very  nearly 
as  high  as  the  covariance  scores.  None  of  the 
techniques  Improved  substantially  upon  the 
scores  of  persistence. 

The  failure  to  Improve  over 
persistence  casts  doubts  upon  the  suitability 
of  simple  extrapolation  as  a  short-range 
forecast  technique.  An  analysis  was  made  of 
cloud  structure  and  evolution  at  several 
locations.  Several  complicating  factors 
were  found,  including  moving  and  stationary 
patterns  in  close  proximity,  sometimes  at 
the  same  level  in  the  vertical,  and  sometimes 
at  different  levels. 

While  it  waa  disheartening 
to  see  that  simple  extrapolation  waa 
acarecely  better  than  persistence,  the 
probelms  can  be  diagnosed  with  the  aid  of 
the  geosynchronous  satellite  data.  Since 
theprablema  can  be  diagnosed,  che  construction 
of  a  successful  short-range  forecast  procedure 
should  be  possible. 
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